Binary spreading process with parity conservation 



Kwangho Park 1 ' 3 , Haye Hinrichsen 2 , and In-mook Kim 3 
1 Theoretische Physik, Fachbereich 10, Gerhard-Mercator-Universitat Duisburg, 47048 Duisburg, Germany 
2 Theoretische Physik, Fachbereich 8, Universitat GH Wuppertal, 42091 Wuppertal, Germany 
3 Department of Physics, Korea University, Seoul, 136-701, Korea 
(February 1, 2008) 

Recently there has been a debate concerning the universal properties of the phase transition in 
the pair contact process with diffusion (PCPD) 2A —* 3A, 2A — > 0. Although some of the critical 
exponents seem to coincide with those of the so-called parity-conserving universality class, it was 
suggested that the PCPD might represent an independent class of phase transitions. This point 
of view is motivated by the argument that the PCPD does not conserve parity of the particle 
number. In the present work we pose the question what happens if the parity conservation law is 
restored. To this end we consider the the reaction-diffusion process 2A — » 4A, 2A —* 0. Surprisingly 
this process displays the same type of critical behavior, leading to the conclusion that the most 
important characteristics of the PCPD is the use of binary reactions for spreading, regardless of 
whether parity is conserved or not. 



PACS numbers: 05.70.Ln, 64.60.Ak, 64.60. Ht 

In the field of nonequilibrium critical phenomena, the 
study of phase transitions from fluctuating into absorbing 
states continues to attract considerable attention jl| . It is 
believed that phase transitions into absorbing states can 
be categorized into a finite number of universality classes 
characterizing the long-range properties at the critical 
point. 

So far two universality classes are firmly established. 
The first and most prominent one is the universality class 
of directed percolation (DP) [§J|], which describes the 
spreading of particles according to the reaction diffusion 
scheme 

A^2A, AA0, (1) 

where A and /x are the rates for offspring production and 
particle decay, respectively. In addition, particles are al- 
lowed to diffuse and the maximal density of particles is 
limited. Therefore, if A is sufficiently high, the system is 
in a fluctuating (active) high-density phase, while for low 
values of A it reaches the (inactive) vacuum state within 
exponentially short time. 

The second established universality class is the so- 
called parity-conserving (PC) class of phase transi- 
tions [^]^6| which appear in spreading processes with 
parity-conserving dynamics such as 

A^3A, 2A^%. (2) 

In this type of spreading process particles can only an- 
nihilate in pairs so that the absorbing phase is charac- 
terized by an algebraic decay of the particle density with 
time. In one spatial dimension, parity conservation al- 
lows the particles to be considered as kinks between op- 
positely oriented domains [QJ|] . Using this interpretation, 
the process can be regarded as a directed percolation pro- 
cess with two ^-symmetric absorbing states (DP2) [[)). 
To some extent, the situation is similar to the one in 
the kinetic Ising model, although in the present case the 



transition is generated by interfacial noise instead of bulk 
noise ju|. 

Apart from these two established universality classes, 
there are only few other possible candidates, the most 
mysterious being the pair contact process with diffusion 
(PCPD), sometimes also called annihilation- fission pro- 
cess. This process was originally introduced by Howard 
and Tauber as a model interpolating between 'real' and 
'imaginary' noise ]Tl[ | and corresponds to the reaction- 
diffusion scheme 

2A^3A, 2A^%. (3) 

Interestingly, this model exhibits a nontrivial phase tran- 
sition even in one spatial dimension. As in the PC class, 
particles can only annihilate in pairs, so that the par- 
ticle density in the inactive phase decays algebraically. 
Moreover, the model has two absorbing states, namely, 
the empty lattice and the state with a single diffusing 
particle. Because of these similarities and an apparent 
numerical coincidence of certain critical exponents, Car- 
Ion et al. raised the possibility that the transition in the 
PCPD might belong to the PC universality class Q. A 
different point of view was presented in Ref . jl3| , suggest- 
ing that the broken parity conservation law in the reac- 
tion diffusion scheme (||) should drive the system away 
from the PC class, leading to the conjecture that the 
transition in the PCPD might belong to a novel, yet un- 
explored universality class. 

Subsequent high-precision simulations ]T^| confirmed 
that some of the critical exponents, especially the or- 
der parameter exponent (3, seem to be incompatible with 
the PC hypothesis, supporting the viewpoint of Ref. [fl3|| . 
On the other hand, the simulations revealed unexpected 
difficulties, in particular unusually strong corrections to 
scaling. It turned out that even after 10 7 time steps it 
is not yet clear whether the 'true' scaling regime has al- 
ready been reached. Therefore, the critical exponents 
in one spatial dimension could only be determined with 
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PC 


0.92(2) 


1.83(3) 


3.22(6) 


PCPD 


< 0.6 
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cyclic model 


0.38(6) 


1.0(1) 


1.8(1) 


present work 


0.50(5) 


1.17(7) 


2.1(1) 



TABLE I. Estimates of the critical exponents for directed 
percolation, the parity-conserving class, and various binary 
spreading processes. 



considerable uncertainty (see Table I). 

Concerning the PCPD transition, there are many open 
questions: Do the critical properties depend on the de- 
tails of the dynamics or do they indeed represent an in- 
dependent universality class which has not been inves- 
tigated before? Does the simple scaling picture, which 
involves only a single length scale, still apply or is it nec- 
essary to consider the possibility of multiscaling? What 
is the origin of the scaling corrections and what are the 
precise values of the critical exponents? 

A possible phenomcnological explanation of the transi- 
tion in the PCPD was proposed in This explanation 
is based on the assumption that the most salient features 
of the process are the interplay of (a) diffusing solitary 
particles, and (b) spreading when at least two particles 
meet at neighboring sites. It was conjectured that a cycli- 
cally coupled model with two particle species consisting 
of a DP process and an annihilation process should dis- 
play the same critical behavior as the PCPD. In fact, 
numerical estimates of the critical exponents seem to be 
compatible with the PCPD results. However, the mere 
numerical coincidence within rather large error bars can- 
not yet be regarded as a proof. 

Since in Ref. |l3| the main argument against the PC 
hypothesis has been the broken parity conservation law 
in the PCPD, it would be interesting to find out what 
happens if the conservation law is restored. This can 
be done by modifying the particle creation process in 
the reaction-diffusion scheme (||), e.g., by considering the 
process 



2 A -> AA , 



2A 



(4) 



In this process, the number of particles is conserved mod- 
ulo 2. As a surprising result, which will be presented be- 
low, we find that this modification does not change the 
type of critical behavior at the transition, i.e., the process 
still behaves in the same way as the PCPD, as already 
observed in the corresponding bosonic field theory [ jTl| . 
Thus, in the attempt to understand the physics of the 
PCPD, it would be misleading to focus exclusively on 
the parity conservation law, rather it is more important 
whether we are dealing with a unary or a binary spread- 
ing process. In a unary spreading process (e.g. in DP 
and PC models), a single particle is able to produce one 
or several offspring. Contrarily, in a binary spreading 
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FIG. 1. The density of particles, p(t), times t 236 as a 
function of time for p = 0.0893, 0.0894, 0.0895, 0.0896, and 
0.0897 from top to bottom, averaged over 2000 runs on a 
system with 2048 sizes. 

process such as the PCPD, two particles are required to 
meet at the same or neighboring sites in order to generate 
offspring. 

A. Definition of the model: The process defined 
in (^), which will be studied in the present work, is a 
binary spreading process with parity-conserving dynam- 
ics. It is defined on a one-dimensional lattice with L sites 
and periodic boundary conditions, where local variables 
Si = 0, 1 indicate whether the site is empty or occupied by 
a particle. The model is controlled by a single parameter 
p and evolves by random-sequential updates according 
to the following dynamic rules. For each update a site i 
is randomly selected and a random number z 6 (0, 1) is 
drawn from a flat distribution. Then the following moves 
are carried out: 

- If p < z and site I is occupied by a particle, it hops 
randomly to the left or to the right. If the selected 
target site is already occupied, both particles anni- 
hilate instantaneously. 

- If p > z and the two sites i and i + 1 are occupied, 
this pair of particles generates two offspring to the 
left (sites i— 2,i — 1) or to the right (sites i + 2, i + 3) 
with equal probability. If the generated particles 
land on an already occupied site, they annihilate 
instantaneously. 

As usual, L update attempts correspond to a time incre- 
ment of 1. The dynamic rules given above can also be 
defined in terms of a reaction-diffusion scheme 
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FIG. 2. Finite size data collapse according to (6) for sys- 
tem sizes L = 64, 90, 128, 180, and 256 averaged over 50000 
runs. 

where q = 1 — p. 

B. Numerical analysis: In order to estimate the crit- 
ical exponents characterizing the transition between the 
active and the absorbing phase we performed standard 
Monte Carlo simulations. To this end we measured the 
density of particles, p(t) — j^^i Si(t), starting with a 
fully occupied lattice as initial condition. At the crit- 
ical point, this quantity should decay algebraically as 
p(t) ~ t~ s . Using this criterion (see Fig. [j]) we estimated 
the critical point by p c — 0.0895(2). For the decay expo- 
nent we obtain the estimate 

5 = p/ V \\ = 0.236(10). (5) 

Next, in order to determine the dynamic exponent z — 
v\\/vx, we performed finite size simulations at the criti- 
cal point. Here the density of particles should obey the 
following finite-size scaling form 

p(t,L)~t- s f(t/L z ), (6) 

where / is a universal scaling function. Using the previ- 
ous estimate 5 = 0.236, the best collapse is obtained for 
z = 1.80(5) (see Fig. ||). Similarly, the third independent 
exponent vh can be determined by studying the behav- 
ior of the density of particles below and above criticality. 
Here we expect the scaling form 

p(t,e)~t- 6 g(t&), (7) 

where e = \p — p c \ denotes the distance from the critical 
point. Using the estimate 8 = 0.236, the best collapse is 
obtained for v\\ = 2.1(1) (see Fig. ||). Combining these 
estimates we arrive at the result 

P = 0.50(5), v± = 1.17(7) , i/M=2.1(l). (8) 
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FIG. 3. Data collapse for off-critical simulations according 
to the scaling form (7) for e = 0.0001, 0.0002, . . . , 0.0064 
averaged over 2000 runs. 

As an additional test, we performed dynamic simulations 
starting with a seed of a single pair of particles located in 
the center, measuring the survival probability P(t) that 
the system has not yet reached an absorbing states, the 
average number of particles N(t), and the mean square 
spreading from the origin R 2 (t) averaged over the sur- 
vival runs. At criticality these quantities should obey 
the power laws P(t) ~ t~ 5 ' , N(t) ~ f, and R 2 {t) ~ t 2 /\ 
where 5' and r\ are dynamical exponents. Measuring 
these quantities we observe strong corrections to scaling, 
which make it impossible to obtain reliable estimates for 
the critical exponents. Fitting straight lines over the last 
decade we find the values (see Fig. Q) 

5'wO.l, ijwO.2, 2/2 w 1.15 (9) 

without being able to estimate the errors. Nevertheless 
the estimate for 2/z is in rough agreement with the pre- 
vious estimate z = 1.80(2). 

We also measured the density of pairs of particles, 
which can be used as an alternative order parameter in 
the present model. Here we find the same type of critical 
behavior, although with slightly different estimates for 
the critical exponents. 

C. Discussion: As shown in Table 1, our estimates 
for the critical exponents are in fair agreement with those 
of the PCPD. In particular, we can rule out the possibil- 
ity of a DP or a PC transition. This result is surprising 
since it suggest that we can introduce an additional sym- 
metry without changing the critical behavior of the tran- 
sition. This means that parity conservation is irrelevant 
for the long-range properties at the transition. 

To understand this observation, we note that there is 
another well-known example where parity conservation 
is irrelevant, namely, the annihilation process 2A in 
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FIG. 4. The survival probability P(t), the average num- 
ber of particles N(t), and the mean square spreading R 2 (t) 
starting with a single pair of particles. 



comparison to the coagulation process 2A — ► A. Both 
processes are known to belong to the same universality 
class and can even be related by an exact similarity trans- 
formation This is due to the fact that the even and 
the odd sector in the parity-conserving process 2A — » 
are essentially equivalent since in both cases the particle 
density decays algebraically until the system is trapped 
in an absorbing state (namely, the empty lattice or a 
state with a single diffusing particle). Breaking the par- 
ity conservation law by a weak perturbation the system 
begins to switch between the even and the odd sector. 
However, this ongoing switching process does not change 
the universal behavior since - from a macroscopic point 
of view - the physical properties of both sectors can- 
not be distinguished. In the present model the situation 
is quite similar. In both sectors we have a transition 
from an active phase into an absorbing state. Therefore, 
the physical properties of both sectors are essentially the 
same so that the breakdown of parity conservation does 
not change the critical behavior. 

In the PC class, however, the situation is completely 
different. In this case parity conservation is indeed rele- 
vant. For example, in the branching-annihilating random 
walk with even number of offspring A — * 3A, 2A — > the 
two sectors are not equivalent because only one of them 
has an absorbing state. Therefore, even a tiny violation 
of the conservation law drives the transition away from 
the PC class. 

How can we verify whether parity conservation in a 
given system is relevant or not? One way would be to 
investigate how the critical behavior changes if the sym- 
metry is broken. Another much more elegant method 
would be to compare seed simulations in the even and 
the odd sector (i.e., starting with two or three particles). 
Here the survival probability P(t) has to be defined as 
the probability that there are at least two particles left. 
If the survival exponent 5' and the critical initial slip 
exponent r\ are different in both sectors (as they are in 
the case of the PC class), parity conservation is relevant. 
However, if the exponents do not depend on the sector 
(as in the present model), we expect the parity symmetry 
to be irrelevant. 



In the light of these results, the conjecture of Ref. E3] 
has to be refined. It is true that we cannot have PC 
critical behavior in systems without parity conservation 
or an equivalent Zi symmetry. On the other hand, the 
broken parity conservation law is not the main charac- 
teristic of the transition in PCPD, rather it is possible 
to restore this symmetry without changing the critical 
behavior. Therefore, a necessary condition for existence 
of this class, for which our understanding is still incom- 
plete, seems to be the binary nature of the process for 
offspring production, i.e., two particles have to meet at 
the same place in order to create new particles. 
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